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Abstract 

Cross sections resulting from scattering that proceeds via an intermediate resonance are shown to 
be exceptionally controllable using a coherent superposition of only two initial states. Full quantum 
computations on F + HD(d = 0; j = 0, 1) — > H + DF,D + HF, which exhibits a resonance in one 
of the reactive channels, support the formal arguments, showing that control is indeed vast. In 
this case the ratio of reactive integral cross sections can be altered by a factor of 62 (compared 
to a non-coherent factor of only 3.3), while the ratio of reactive differential cross sections can be 
altered by a factor of over 6000 (compared to a non-coherent factor of less than 7). These results 
constitute the first prediction of extensive quantum control in a collisional process. 

PACS numbers: 34.50.Lf, 34.50.Rk 
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Coherent control is a rapidly developing method for controlling atomic and molecular pro- 
cesses. Recent experimental and theoretical developments in this area have been extensively 
summarized in Refs. [1, 2]. In this approach one manipulates, through laboratory parame- 
ters, quantum interference contributions that arise when a final state is reached by two or 
more coherent routes. Although initially introduced for unimolecular processes, coherent 
control has been extended formally to collisions [3-6] and has been numerically studied for 
reactive scattering of the isotopic variants of H + H 2 [3] and for H 2 + H 2 [6] . All of the 
established principles apply equally well to atomic, molecular and nuclear processes. 

There are very few general results known for control of scattering. We previously showed 
that complete control, i.e. the ability to completely suppress or maximally enhance a par- 
ticular product arrangement channel, is achieved in collision processes when the number 
of reactant states used in an initial coherent superposition is greater than the number of 
open final product states in that product channel [5]. In this article we show that final 
product channels that arise via the formation of a resonance can be completely controlled 
via a superposition of only two initial states, regardless of the number of open product states. 
We provide, as a numerical example, results of a full 3D quantum computation of control 
through the recently discovered reactive resonance in F + HD — > D + HF[7-9]. Even though 
this resonance is short-lived [8], and direct reactive scattering occurs as well, we find that a 
very large range of control is possible. 

Control of a collisional process is achieved[3] by constructing an initial state |\&) in the ini- 
tial arrangement channel composed of a superposition of N energetically degenerate asymp- 
totic states \q,m). Here q denotes the arrangement channel and m indicates all state labels; 
the energy label is dropped for convenience. For the case of N = 2: 



The cross section a q > for scattering into final arrangement channel q', starting from state 
is given by 



where T is the transition operator. For notational convenience we drop the q 1 subscript 
unless necessary. 

Equation (2) can be rewritten as 



|tf> = Cl |g,l) + c 2 |g,2) . 



(1) 
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where 

Mu) = Y,^AT\q\n){q'MT\q, 3 ) . (4) 

n 

To produce the integral cross section, the sum over n includes an integral over the product 
scattering angles as well as a sum over all final diatom states. 

Equation (3) can be rewritten by defining the relative amplitude s = 
| c 2 1 2 / [ I Ci 1 2 + | c 2 1 2 ] , < s < 1 , and the relative phase between the two coefficients as 
012 = Arg [c 2 /ci]. The cross section then becomes 

oy = (1 - s)oy(ll) + sa q ,(22) 



+ 2y/s(l-s) |<v(12)| cos (Arg[^(12)]+0 12 ) , (5) 

indicating that control is achieved experimentally by varying s and 0i 2 . Note that ay (11) 
and ay (22) are real, corresponding to the cross sections for initial states 1 and 2 respec- 
tively, while a q i{12) is complex, corresponding to quantum interference between the two 
reaction pathways. Since 12 can always be chosen so that Arg[(T g /(12)] + 12 is a multiple 
of 2-7T, the extent of control over a q i is determined by |ay(12)|, which is seen to satisfy the 
Schwartz inequality |a g /(12)| < ^ a q i{ll)a q i{22). Maximum control therefore occurs when 



|ay(12)| = yj a q '(ll)a q/ (22). The value of the cross section can be controlled within the range 
<T g /(min) < o q i < (Tg/(max), where the limiting values and corresponding control parameters 
are [assuming a g /(22) > o- g /(ll)] 

(Tg/(min) = ff,'(max) = a q i{\\) + a q i{22) 
s mm = — 7-n- t^t s max — 2 



v (ii)+ v (22) V (ii)+ V (22) 
0i 2 (min) = 7T - Arg[oy(12)] 0i 2 (max) = -Arg [^(12)] . 

Note that <r g /(max) is twice that achievable in the absence of the interference term. 

We show below that the Schwartz equality is satisfied for scattering into a particular 
arrangement channel if it occurs through a resonance at all scattering angles and for all 
product states, a situation that we term an ideal resonance. By comparison, the limits 
of non-coherent control (i.e. when no quantum interference occurs) are the cross sections 
corresponding to the two initial states <r g /(ll) and <T g /(22) (corresponding to s = and 

5 = 1). 

Consider a system displaying an isolated Feshbach resonance [10]. In this case the T 
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matrix elements at scattering energy E are of the general form[10, 11] 

(q',n\T\q,j)= 1 ' r (q',n) lr (q,j)/[E-E r } , (6) 

where Y r (q',n) and 7 r (9,j) are matrix elements of the Hamiltonian coupling the product 
state to the resonance state (denoted by the subscript r) and coupling the resonance to the 
initial state, respectively E r is the complex energy associated with the resonance. If the 
form of Eq. (6) holds, for all final states and all scattering angles at the energy of interest, 
we have 



n 



„ / 19 x 7r*(g, l)7r(g,2) s^^u ^M2 



so that 









£/ — E r 


2 



0y = — — =rr^- X 



(1 - s)| 7r (g, 1)| 2 + s| 7r (g, 2)| 2 + 2yJs(l-sM(q, l) 7r (g, 2)| cos (Arg[^(12)] + 12 ) 



(5) 



Under these circumstances the magnitudes of the cr q >(ij) satisfy the Schwartz equality and 
complete control over ay is possible. 

Control via a resonance over the ratio of cross sections into different arrangement channels 
q', q" requires that at least one of these channels have a direct scattering component, or that 
they arise by scattering through different ideal resonances. If this is not the case, i.e. if 
both final arrangements q' and q" are accessible only via the same ideal resonance, then, as 
is evident from Eq. (8), the ratio a q > ja q n no longer depends upon the control parameters s 
and 0i2- 

The above argument is general and the recently observed reactive resonance in F + HD — > 
D + HF scattering [7-9] provides an opportunity for a challenging numerical test within 
the framework of molecular scattering. The FHD resonance, corresponding to a collinear 
arrangement with quantum numbers (vpn = 3,t>HD = ®,Vbend = 0), is short-lived, with a 
lifetime of 109 fs. In the low energy regime, where the resonance dominates the reactive 
scattering, the rotational period is 1420 fs, and the contributions from several total angular 
momentum values overlap one another [8]). Advantageously, the collinear FHD complex is 
much more likely to decay to D + HF than to H + DF, so that complete control over the 



resonant mechanism should allow large control over both the D + HF cross sections as well 
as over the ratio of D + HF to H + DF reactive cross sections. Deviations of control from 
the maximum value can also be be enlightening, being attributable to non-resonant, direct 
reactive scattering contributions. 

Control results shown below were obtained using a converged full three-dimensional quan- 
tum coupled-channel hyperspherical coordinate approach [12] on the Stark- Werner potential 
energy surface [13], using an available code [14]. Results were obtained for a variety of ener- 
gies over the range 0.25-0.31 eV (relative to the HD minimum) in increments of 0.005 eV. A 
basis set containing all three channels with maximum diatomic energy levels of 1.7 eV, max- 
imum rotational quantum number j max = 15, and maximum total angular momentum and 
helicity J max = 31 and K max = 4 were used. S'-matrix elements for all possible combinations 
of energetically available initial and final states were calculated. 

Control was studied for the case where the initial state was comprised of a superposition 
of the (t>i,ji,mi) = (0,0,0) and (t>2, J2, ^2) = (0,1,0) states of HD, which have internal 
energies of e\ = 0.23252 and e 2 = 0.24358 eV, and with k\ and k 2 satisfying the conditions 
[4] 

K 1 = K 2 ; E = h 2 k 2 1 /2fi + e 1 = h 2 kl/2fx + e 2 . (9) 

Here ftKj and hki are the center of mass and relative momenta respectively, and \i = 
M F M HD /M FilD is the reduced mass. For the case we have considered the kinetic energies 
of the two coherent components, defined as E kj = h 2 k 2 /2fi are related by E k2 = E kl — 
0.01106 eV. 

Most relevant is the ratio of product cross sections. Figure la shows the maximum and 
minimum values of the ratio of integral reactive cross sections, r = a[D + HF]/er[H + DF] 
attainable at each energy. Results using coherent control are shown as solid lines while those 
without control [i.e. resulting from ay (11) and ay (22)] are shown as dashed lines. Clearly, 
r peaks near the resonance (E = 0.2550 eV) with the ratio controllable, at resonance, over 
a huge range: 2.60 to 161, a factor of 62. This range is far greater than the analogous 
non-coherent control factor of 3.3. Further, the coherent control parameters resulting in the 
minimum ratio (2.60) are s = 0.7027 and </> 12 = 185.1°, while for the maximum (161) they 
are s = 0.9344 and 0i 2 = 352.9°, sufficiently far apart to allow easy discrimination. 

These results are shown in a somewhat different way in Figure lb, where the ratio of 
maximum and minimum ratios r, i.e. R = max[r]/ min[r] attainable via coherent control 
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are compared to those attainable without coherent control. Three features are evident: 
coherent control affords a vastly larger range of control, the best ratio control is seen at the 
resonance, and R for the coherent control case remains greater that 20 at energies up to 
0.31 eV, where the cross section is no longer dominated by the resonant process. However, 
as seen in Fig la, in this higher energy region the actual value of r is quite small. 

Figures lc and Id provide the cross sections from which the earlier panels are constructed. 
Figure lc shows that the D + HF cross section can be coherently controlled at the resonance 
from a minimum of 0.0850 A to a maximum of 2.193 A. This is to be compared to the 
ideal case of complete control with a minimum of and a maximum, corresponding to 
ay(ll) + ay(22), of 2.278 A. Alternatively, we can gauge the "quality" of the resonance from 
the deviation of ay(12) /[a q i(ll)a q i(22)] 1 ^ 2 , from unity. This ratio peaks, at the resonance 
energy, at a value of 0.906. As a general observation we note that this ratio is far larger 
than that seen in previous studies, such as those of the isotopic analogs of H + H 2 [4]. Hence 
control is far more extensive in this case than in any collisional process previously studied. 

Results were also obtained for control of the differential cross section. In this case the 
sum in Eq. (4) does not contain an integral over scattering angle. Results at the resonant 
energy show that the greatest control occurs for backward scattering (8 = 180°). There the 
ratio of reactive cross sections can be coherently controlled from a minimum of 0.236 (at 
s = 0.6469 and <p l2 = 167.0°) to a maximum of 1455 (at s = 0.8097 and <p l2 = 331.1°), a 
factor of over 6000. By contrast, the non-coherent control factor never rises above 7 at this 
energy. At this scattering angle the Schwartz ratio a q '(12) /[a q i(ll)a q i(22)] 1 ^ 2 = 0.984. Such 
near-unity values are consistent with the fact that much of the averaging is eliminated in 
examining the differential cross section. 

We have shown above that if there is an ideal resonance then complete control is possible 
using a superposition of two initial states. We note further that the behavior of the Schwartz 
ratio at the resonance energy provides a useful theoretical test for the contrapositive. That 
is, if the Schwartz equality does not hold at a particular energy then the dynamics must 
have a non- resonant component. Thus, for example, information regarding the resonant 
(or non-resonant) behavior of H + HD scattering in the forward direction, a subject of 
recent debate[15], can be augmented by computing the Schwartz ratio in that direction. 
Doing so using the same initial superposition as for the F + HD above and the established 
BKMP2 potential surface[16] gives us a q/ (12)/[a q/ (ll)a ql (22)] 1 ^ 2 = 0.850, indicating that the 
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dynamics does not go solely through a resonance. 

In conclusion, we have shown that scattering through a resonance allows for a vast range 
of control over product cross sections using an initial state comprised of a superposition of 
as few as two states. Preparing this superposition for atom-diatom scattering does present 
a challenge since the terms in the superposition are comprised of wavefunctions where the 
translational and the internal diatom wavefunctions are correlated [4]. However, significant 
advances have been made in experimental coherent control studies for unimolecular prob- 
lems [1, 2] and there are a number of possible approaches to preparing the states required 
for bimolecular control. These include extending a method successfully applied [17] to atoms 
in order to prepare the required correlated initial states, or state-selectively accelerating or 
decelerating molecules[18] that are in different internal state to achieve this result. Alter- 
natively, we note that, for the case of the scattering of identical particles, preparation of 
the required states is straightforward [6]. Since the discussion in this paper applies equally 
well to such cases we suggest that the effect described here may well be first observed, in 
molecular physics, in the scattering of identical diatomic molecules, e.g. OH + OH. 

We thank Dr. D. Skouteris for sending us his auxiliary code for calculating cross sections 
from the scattering output. This work was supported in part by the U.S. Office of Naval 
Research and by the Natural Sciences and Engineering Research Council of Canada. 
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FIG. 1: Maximum and minimum values through which integral cross sections can be coherently 
(solid line) or non-coherently (dashed line) controlled using the initial state described in the text, 
(a) The ratio of the cross sections into the D + HF channel to that in the H + DF channel. Note 
that the minimum values are very close to the r = line; (b) The ratio of the maximum to 
minimum r at each energy; (c) the cross section for the resonant D + HF channel; (d) the cross 
section or the non-resonant H + DF channel. 
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This figure "fig.gif" is available in "gif" format from: 



http://arXiv.org/ps/quant-ph/0405010vl 



